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In this paper, we propose a new class of parametrization of the equation of state of dark energy. 
In contrast with the famous CPL parametrization, these new parametrization of the equation of 
state does not divergent during the evolution of the Universe even in the future. Also, we per- 
form a observational constraint on two simplest dark energy models belonging to this new class 
of parametrization, by using the Markov Chain Monte Carlo (MCMC) method and the combined 
latest observational data from the type la supernova compilations including Union2(557), cosmic 
CNj . microwave background, and baryon acoustic oscillation. 
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I. INTRODUCTION 

The cosmic acceleration has been discovered since 1998 [l|, HJ, but we are still in the dark about the nature of 
this mystery. The cosmological constant proposed by Einstein to build a static Universe solution is the simplest 
dark energy model to explain the acceleration, but it suffers from the fine-tuning and coincidence problems. On the 
other side, current observations do not exclude definitely the time dependent evolution of dark energy, then many 
alternative dark energy models were proposed in the literatures, for a recent review, see @, 0j- 

The evolution of dark energy is determined by its equation of state (EoS), which can be obtained in a particular 



model or constructed directly from observational data. In the later case, we usually make some assumptions on the 
form of the EoS as a parametrization for the dark energy model, then extract the information from observations. A 



most wildly used parametrization is called the Chevallier-Polarski-Linder (CPL) model [5|, |6[ with the form of EoS as 

w(z) = w + wi— — , (1) 

>'. 

with redshift z. Here, wq is the present value of the EoS and w\ is the derivative with respect to the scale factor. 
\Q • Comparing with the linear form w(z) = wo + w\z, the CPL parametrization is well behaved and bounded at high 
redshifts. This is the direct motivation of proposing such a form (flj. However, the CPL form also suffers from the 
divergence problem when z approaches —1, which corresponds to future evolution of the EoS. This is undoubtedly 
\Q 1 a nonphysical feature, and thus it prevents the CPL parametrization from genuinely covering the scalar field models 
as well as other theoretical models. Recently, a logarithm form w(z) = wq + w\ Jj_n(2 + z)(l + — In 2] and an 
oscillating form w(z) = wo + w-y [sin(l + z)(l + z)^ 1 — sin l] are proposed in Ref. to overcome the shortcoming 
of the CPL model. 

To extend parametrization of dark energy to redshifts z — > — 1, we suggest two new models 



X 



Model I: w(z) = w + wi ^ , (2) 

1 + z l 

and 

z 2 

Model II: w(z) — wq + w%- ^ , (3) 

which arc free of divergence in the whole range of the redshift z 6 [— 1, oo). Clearly, Model I reduce to the linear 
form when z < 1 as the CPL form, while Model II will be not. Actually, one can even generalize the form to a more 
general case while keeps the advantage of divergence-free, as long as some conditions are satisfied. We will make some 
discussions about it in the conclusion section, and one shall see that the forms @ and ([3]) are the most simplest ones 
of this class of parametrization. 
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We perform a global data fitting analysis on these two divergence-free parametrization of dark energy models and 
constrain the parameters from the current observational data, including the full CMB power spectrum from seven- 
year WMAP data (WMAP7) [13], the baryon acoustic oscillations (BAO) data from SDSS DR7 p~D] and Union2 
SNIa with 557 samples [13]. When doing the analysis, we also take the perturbation of dark energy into account, 
because the perturbation plays a important role when one confronts the dynamical dark energy models with the 
current observations [13H16I ]. For instance, when the dark energy dominates the Universe at later time in a spatial flat 
model, the gravitational potential due to the perturbation of dark energy becomes time-dependent, thus there will be 
a significant imprint on the large-scale CMB power spectra by the integrated Sachs- Wolfe (ISW) effect [17H22| . We 
also make some analysis on this effect in this paper. We added two new parameters Wo,wi and modified the public 
available CosmoMC package [23[ to satisfy these two models. During the calculation, we run 32 independent chains 
for each case, and we make sure the convergence of the chains by typically getting 1Z — 1 to be less than 0.01, see 
Ref. [23M27I ] and the readme file on the website of CosmoMC. 

This paper is organized as follows: In SecJITl we will study the background and perturbation evolutions for both 
models. In Sec lIII| we perform a global fitting to the current observational data by using the Markov Chain Monte 
Carlo (MCMC) method. Some analysis on the ISW effect of there two models are given in Sec lIVI In the last 
section, we will give some conclusions and make some discussions on how to generalize these models without losing 
the advantage of divergence-free. 



II. BACKGROUND AND PERTURBATION EVOLUTION 



As we mentioned before, the new parametrizations of dark energy ([2|) and ([3]) are free of divergence in the whole 
range of the redshift z g [—Loo), especially when z — » —1. For clearness, we list the limiting values of w(z) when 
z -> 00, -1, in TabUfor the CPL model, Model I and Model II. At low shifts, the Model I form will reduce to the 
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TABLE I: The limiting values of the EoS of CPL, Model I and Model II. 



linear form w(z) wq+wiz as the CPL model, while Model II will reduce to the quadratic form, i.e. w(z) « wo+wiz 2 , 
and both of them are free of divergence. From the continuity equation of the dark energy with EoS (O and one 
can obtain the its energy density, and then the Friedmann equation would be 
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fljkf (1 + z? + ^r(l + zf + (1 - n 
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(4) 



for a spacial flat Universe. Here we have defined fii = pi/p c , {i = b,m,r) and Qm = + ^m- The functions f(z)± 
are given by 



f±{z) = exp 



± — — arctan(z) 



(l + z)T— (1 + z 



2\ — 



(5) 



where / + and /_ correspond to Model I and II, respectively. 

As we known, for dynamical dark energy models, namely, the EoS of dark energy is not constantly equal to —1, it is 
crucial to include dark energy perturbations [13l4l6| . Working in the conformal Newtonian gauge with the perturbed 
metric as the following 



ds 2 = a(?7) 2 - (1 + 2iP)dr] 2 + (1 - 2cj))S ij dx i dx j 

one can obtain the perturbed equations of dark energy as [2814301 ] 

5' = -(l + w)(0 -i<p')-m{c 2 s -w)5 , 



-H(l - 3w)6- 



9 -\ — 

1 + w 1 + w 



k 2 S + k 2 ^, 



(6) 

(7) 
(8) 
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where prime denotes the derivatives with respect to conformal time 77, H = a' /a is the conformal Hubble parameter. 
Here 5 and 9 are the density and velocity perturbations respectively. The square sound speed c 2 is defined as 
c 2 = Sp/Sp. It is clear that for the quintessence-like (or phantom-like) models with w > — 1 (or w < —1), the 
perturbation equations ([7} and © are well-defined. In other words, as long as the EoS does not cross the cosmological 
boundary, the perturbation equations will not divergence. However, when the cross happens in the dark energy 
models such as CPL, Model I and II, one should take care when w — — 1. This is so-called the divergence problem 
for perturbations of dark energy. A very useful practical solution to this problem is to introduce a small positive 
constant e, then the whole range of the allowed values of EoS is divided into three parts: w > — 1 + e, w < — 1 — e 
and — 1 — e <w < — 1 + e, see Ref. [l3j]. For the first two regions, the perturbation equations are well-defined without 
divergence. To avoid the divergence in the last region, we take the approach that match the perturbations in this 
region to the first two regions at the boundary and set 5' = 9' = in this region. In our numerical calculations, we 
limit the range to be e < 1CP 7 , suggested in Ref. 311 ] . 



III. GLOBAL FITTING PROCEDURE AND RESULTS 

In this section, we will use the Markov Chain Monte Carlo (MCMC) method to make a global fitting on the 
cosmological parameters in Model I and II. We have modified the public available CosmoMC package [23:] to satisfy 
Model I and II, in which we have added two new parameters Wq and W\ and the perturbation equations for dark 
energy. Our most general parameter space vector is 

,2 



V= |fi 6 ^,r! c ^,e,T,u;o, W i,n s ,log[10 lu A 5 .]|, (9) 

where ^h 2 and il c h 2 are the physical baryon and cold dark matter densities, O is the ratio (multiplied by 100) of 
the sound horizon to the angular diameter distance at decoupling, r is the optical depth to re-ionization, wq and w\ 
are the parameters of dark energy EoS given by Eqs.® and ([3]), n s is the scalar spectral index and A s is defined 
as the amplitude of the initial power spectrum. The pivot scale of the initial scalar power spectrum we used here is 
k sa = 0.05Mpc _1 and we have assumed purely adiabatic initial conditions. In this paper, we will use the WMAP 
seven-year data [13] , the BAO data [TTj] and also the current supernova " Union2" data with 557 simples [l2[ • All the 
calculations are performed on the Astrophysical Beowulf Cluster (AstroBC), which is a parallel computing network 
system designed, set up and tested by ourselves, and it is much fast and stable than a single personal computer. 

In order to determine the best value of parameters (with la error at least ) in Model I and II, we will use the 
maximum likelihood method and take the total likelihood function C = e~ x I 2 as the products of these separate 
likelihood functions of each data set, thus we get 

X 2 = X 2 CMB + X 2 BAO + X 2 SN ■ ( 10 ) 

Then, one can get the best fitting values of parameters by minimizing \ 2 - 

Our global fitting results are presented in TabJlTl in which we list the best fit values with la, 2a errors for each 
parameter in the parameter space vector ([9]), as well as the constraint for the age of the Universe. For comparison, 
we also present the fitting results without dark energy perturbations in Tab ITU and we find that the results with 
perturbations show a better fit to that without perturbations by comparing the minimal values of \ 2 ■ In our analysis, 
we have fixed the sound speed of dark energy in the rest frame c 2 , which related to c 2 in a general frame as [l8l . |30| 

c 2 s = ~c 2 s+ m(l + w)(b 2 s -c 2 a )^5- 1 . (11) 

To derive the above relation, we have used the useful relation between the gauge invariant, rest frame density pertur- 
bation S and the density and velocity perturbations in a general frame, S and 9 

6 = 5 + 3H(l + w)-^, (12) 

and the definition of the adiabatic sound speed 

2 _ P' W ' /-l Q\ 

p' 3n(l + w) 

see Ref. [l8l . |30| . Since it has been shown that the constraints on the dark energy sound speed in the rest frame c 2 
in dynamical dark energy models are still not accurate enough [32j . then one can treats the dark energy model as 
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a effective scalar- field model and sets c 2 s to be unit Q- The results would not be affected by this treatment if c 2 s 
is considered as a free parameter. In this paper, we have taken the full CMB temperature and polarization power 
spectra instead of just taking the WMAP distance prior, including R, I a and z*, which encoding the background 
distance information and can be used to investigate dark energy models with greatly simplified numerical calculations. 
However, it was found that this distance prior is somewhat cosmological model dependent, thus, one would loss some 
CMB information by using it to fit the parameters (Toj |. For instance, as we known, if there was a time-dependent 
gravitational potential at later times, it will contribute to the ISW effect. Therefore, our results would be more 
reliable than that obtained by using just the distance prior. 
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w/. perturb. 
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l.UJOZ_o. 0072-0. 0094 


r 
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n nsi/i+°' 0526 + ' 0679 

U.UOiy;_o. 0319-0. 0409 


n nsR7+°- 0475 + - 0619 

U.UOO / -0.0379-0.0453 


n ns^Q+ '0 449 +0' 0614 

U.UOOO_ . 0333-0. 0473 


wo 


1 m /lo+0. 5183+0. 6482 
l.U14:O_ .5907-0 .7247 


1 m cq+0. 6181+0. 7664 
l.UlOO_ .5576-0.7103 


1 091 /| +0.3093+0.4287 
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Wl 
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fl 1 1/IC+2. 0636+2. 6046 
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n s 
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lO.OU11_ 3608-0. 4255 
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Xmm 


8000.764 
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TABLE II: The fitting results of the parameters with la, 2a regions in Model I and Model II with and without perturbations 
from the combination of WMAP, BAO and SN. 

From TabHIl one can see that in both models, the best-fit values of the scalar spectral index n s are smaller than 
1. That is to say, the scalar spectrum is "red" tiled. The difference between Model I and II is shown in FigHJ in 
which we have plotted the evolution of the EoS and the adiabatic square sound speed for both models with their 
best- fit vafues of u>o and w±. From FigdJ we find that the best- fit dark energy model is a phantom model, whose 
w(z) is always smaller than —1, and the evolution of w(z) is quite different between these two models. However, the 
difference of the present values of w(z) is smaller than one percent. These results imply that though the dynamical 
dark energy models are mildly favored, the current data cannot distinguish them effectively. 

In Fig. [21 we show the CMB temperature power spectrum for Model I and II with their best-fit parameters listed 
in Tab. [TlJ The spectrum for ACDM is also plotted for comparison. We find that the most significant discrepancies 
between Model I, II and ACDM come from the low Is, or large angular scales (~ 180°/^). The contribution to the 
secondary CMB anisotropy is mainly contributed from the late ISW effect, see the discuss in next section. Both cures 
of Model I and II tend to be plat on large scales. In the Fig|3]and FigJH we also show the one dimensional probability 
distribution of each parameters and also the two dimensional contour plots between each other in the Model I and II 
from combination of WMAP, SN and BAO. 



IV. THE ISW EFFECT 



The ISW effect is due to the existence of time-dependent gravitational potentials along the line of sight from the 
epoch of last scattering up to now [13, [H, H3] . There is a net change in the energy of a CMB photon as it passes 
through these evolving gravitational potential wells, and the sum of these contributions is called the ISW effect. The 
ISW contribution to the CMB power in the Ith multipole is given by 

2l + 1 c isw = J_ f dk ki \e[ sw { VQ ,k)\^ 



Ait 1 2tt 2 J k 21 + 1 

where Qj sw is the multipole decomposition of the fractional fluctuation due to the ISW effect 



21 + 1 



= 2 ^'(»?,fc)ji[fcfa)-»7)]<&7. (15) 
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FIG. 1: The evolution of EoS w(z) and the adiabatic square sound speed in Model I and II with their best-fit values of wo 
and wi . 



Here, -q and 77* denote the present and the last scattering conformal times respectively, and ji are the spherical Bessel 
functions. As we known, the perturbation equation of ip = 4ttG(p + p) x / 2 u is given by 



u" - c 2 V 2 u - —u = . 



(16) 



where 8 = (1 + w) 1 l' 2 a 1 . The adiabatic square sound speed in a dark energy dominated Universe is given by 
Eq. (|13[) . and one can rewritten c 2 in terms of redshift 



9 1 + z dw 

c s =w + 



3(1 + w) dz 

where w takes the form of @ and ([3]). And thus, we obtain 



Model I: 



ZWi 



c a (z) = w Q + + 



(l-z)(l + z) 2 Wl 



and 



Model II: c 2 a (z) = w + 



l + z 2 3(1 + z 2 ){\ + z 2 + (1 + z 2 )wq + zwi) 
2z(l + z)wi 



Z 2 Wi 



+ 



1 + z 2 3(1 + z 2 )(l + z 2 + (1 + z 2 )w a + z 2 wi) 



(17) 



(18) 



(19) 



The evolution of c 2 is shown in Fig. [T] in Model I and II. 

It is usual to divide the ISW effect into early and late parts. The early ISW effect corresponds to the variation in the 
perturbation potentials which occurs in the time interval between matter-radiation equality and recombination, while 
the late ISW effect corresponds to the variations caused by the presence of dark energy. The expansion time scale at 
dark energy domination rjd e sets a critical wavelength corresponding to ki]d e = 1. Thus, for adiabatic perturbations, 
the late ISW effect in Eq. (|15|l during the dark energy epoch rj 3> rjde for a given mode can be approximated by [H, H3] 



21 + 1 



(20) 
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FIG. 2: The CMB Cf 1 power spectrum v.s. multipole moment I, where the yellow dots with error bars denote the observational 
data with their corresponding uncertainties from WMAP 7-year results, the red solid line, red dashed and black dotted-dashed 
curves are for Model I, II and ACDM model with their best-fit values from WMAP+BAO+SN joint constraint. 



where rjk = t/q — (I + l/2)/k is the conformal time when a given fc-mode contributes maximally to the angle that this 
scale subtends on the sky, which obtained at the peak of the Bessel function ji . Here the integral // is given by 

m ' 2 r[(Z + 2)/2] V ' 

Substituting the results above into Eq. ([T4"]) . we obtain the Ci momenta of the late ISW effect 

^"'i^MT*''^- (22) 

Therefore, one can at least numerically obtain the solution of ip from Eq. (|16p and then get C/s. From Eq. (|2"2")l , one 
can see that the main contribution of the late ISW effect mainly comes from the low momenta I, see Fig. [2j 



V. DISCUSSION AND CONCLUSION 



In this paper, we have proposed two new parametrizations for the EoS of dark energy, see Eqs.([2]) and ([3]). Actually, 
these two models belong to the following general case 



w(z) = Wq + Wi 



1 + z 



2)1 



(23) 



where Model I corresponds to (to = l,n = 1) and Model II corresponds to (to = 2,n — 1). From the above 
equation, one can see that as long as n is a non-zero integer, w(z) would be not divergent when z — > —1, and the 
limiting value is w(z — > —1) = wo + (— l) m iui/2. Furthermore, if m < 2n, the EoS would be also not divergent 
when the redshift goes to infinite. The limiting values for w(z) when z — and z — > oo are w(0) = wq and 
w(z — > oo) = wo(m < 2n), wq + wi(m = 2n). In general, one can not get an analytic result for the function f{z) in 
the background equation and usually it is given by the following integration 



f(z) = exp 



3wi 



o (l + z^)(l + ~z) 



(24) 
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FIG. 3: The 1-D constraints on each parameters in V = {SI5/1 2 , fl c h 2 , ©, t, wo, tOi, n s , log[10 10 A 3 ] } with the age of Universe, 
and 2-D contours on these parameters with la, 2a confidence level each other in Model I. Here the dotted curves in each f-D 
plots are the mean likelihood of the samples, while the solid curves are the marginalized probabilities for each parameters. 



It is interesting that in Model I, the EoS could be rewritten as 

w{z) = w + - 1 ( 25 ) 

from which one can see that w(z) wq + w\j z when the redshifts are large {z 3> 1), while w(z) w wq + w\z when 
the redshifts are small (z <C 1), in other words, the EoS of Model I has the symmetry of z <H> z^ 1 . 

In conclusion, we have performed a global fit study on two parametrization models for dark energy, named Model 

I and II. In contrast to the famous CPL parametrization, in which w(z) is divergent when z — > —1, Model I and 

II overcome this shortcoming and are divergence- free in the whole range of the redshifts z G [— l,oo). We have not 
only obtained the best- fit values with la and 2a regions for both of Model I and II with and without dark energy 
perturbations, but also found that the differences of CMB temperature power spectrum Cf T on large angular scales 
are significant between Model I, II and ACDM model. There also exists a flat plateau on the large scales in both 
models. In view of this difference, on which the ISW effect plays an important role. The effects on the parametrizations 
of Model I and II made by the ISW data from the correlations of CMB and large scale structure will be studied in 
our further works. And the class of this parametrization (I23[) is also worth investigating in future. 
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